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Abstract: 

I discuss a Monte Carlo algorithm for hadronization in the Lund string fragmentation picture, 
which conserves energy and momentum exactly as is proposed by the model. An embryo 
to this Monte Carlo is used to calculate the total decay density of a string. This is found 
to reach its predicted asymptotic behaviour already at a string energy of about six hadron 
masses. 
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1 Introduction 



In the Lund string fragmentation model the colour field between a quark and an antiquark 
is assumed to be compressed into a 1+1 dimensional string, as the magnetic field in a super 
conductor. The density for a hadronic final state is a phase space factor times an exponential 
suppression of the space-time area swept by the fragmenting string. The model can be 
formulated in terms of an iterative fragmentation process, which makes it very suitable for 
MC implementation. A widely used MC program, Jetset[0], is based on the model. 

As global energy-momentum conservation is ill suited for iterative algorithms, it is observed 
by additional approximations in the MC. In general, these approximations have negligible 
effects, but for subtle enough observables, and investigations at moderate energies, it may 
be worthwhile to point out that there is a difference between predictions of the Lund model 
and the Lund Monte Carlo. 

Due to its success in describing data, the model is being used for more and more detailed 
investigations. One important example is Bose-Einstein correlation investigations at LEP, 
where one appealing model is based on the Lund string area law [||. Another possibility is 
to study of transverse momentum properties in pencil-like two-jet events, which due to the 
large LEPl statistics can be selected without losing statistical significance. Different model 
predictions for these events are discussed in and awaiting confrontation with data. 

In experiments like the B-factories, JLAB and BEPC/BES in Beijing, it is possible to do 
detailed QCD studies at moderate energies. The more moderate is the energy, the more 
important is global momentum conservation treatment in models and Monte Carlo. 

Examples of model vs. Monte Carlo differences, relevant at moderate energy few-body pro- 
duction, are given in ||^, which also presents a method to more precisely generate few-body 
states. A similar approach is used in 0, to calculate (qq— 7r"'"7r~)/(qq— hadrons), for 
moderate-energy qq systems. 

In this paper I present a general algorithm, applicable to any multiplicity, which conserves 
energy and momentum exactly as proposed by the model. A simple embryo of this algorithm 
has been implemented 0, and I use that to present numerical calculations of the total string 
decay density. 

A set of extensions to the algorithm are needed, before a realistic generation of hadronic 
final states are possible. Most of these are postponed to future work, and only briefiy 
discussed in the outlook. The results in this paper are obtained for the simple case of a 
1+1 dimensional string decaying into hadrons of one fiavour only. However, the algorithm is 
presented including two important extensions, namely several fiavours and non-perturbative 
generation of transverse momenta. 

The outline of this paper is as follows: In section |] the Lund string fragmentation model 
is reviewed. Section ^ deals with differences between Monte Carlo and model, and the new 
Monte Carlo algorithm is presented in section ^. The calculation of the string decay density. 
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Figure 1: x — t picture of a fragmenting 1+1 dimensional string. A quark and an anti-quark 
from adjacent vertices form a meson (dotted lines). In the Lund model, the final state is 
exponentially suppressed in the colour coherence area A (shaded region). The vertices tend 
to lie around a hyperbola of constant proper time (dashed line). 

based on this new algorithm, is presented in section ^ and the paper ends with a summary 
and outlook in section |[ 



2 The String Fragmentation Model 

In this review, we start with the simple case of only one hadron flavour and no p_L. The 
density for a string with mass to produce n hadrons is in the model given by 

52(P-^p,)exp(-M), P^ = s. (1) 

r 

Here N and h are free parameters of the model, and A is (the string tension squared times) 
the space-time area swept by the string, cf. Fig 0. The density for a specific n-body state is 
given by the integrand of Qn- The index r runs over the hadrons in rank order., specified so 
that the rank 1 hadron contains the original quark and an anti-quark qi, the rank 2 hadron 
contains the quark qi from the same break-up and an anti-quark q2 from the next vertex, 
etc. 

The total decay density 

9{s) = Y,9n{s) (2) 

n 

can be written as an integral over a production factor for the rank 1 hadron times the decay 
density of the remaining string. For s ^ m? the solution is of the form 

g{s) oc s", (3) 

where a is related to the parameters and hm^ through 

l = N —[l-zYe-^'^''\ (4) 



gn[.s) = 



n „ 

l[N d\6{pl-m^) 

r=l 
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z represents the light-cone momentum fraction of the rank 1 hadron, which thus can be 
picked from a distribution 

d 7 hm ^ 
f{z,m^)dz = N—{l-zyexp{ ). (5) 

z z 

This expression is independent of the string energy, which makes the formahsm suitable for 
an iterative algorithm. Each produced hadron can be seen as the rank 1 hadron of a specific 
string remainder, and its momentum can be determined from Eq (^. Eq (|^) then acts as 
a unitarity constraint. This holds as long as the remaining energy is large enough for g to 
have its asymptotic form in Eq (^), when global energy-momentum conservation effects can 
be neglected. 

The fragmentation function Eq (^) gives left-right symmetry, i.e. we get the same final state 
distributions if we fragment iteratively from the quark end as from the anti-quark end of the 
string. Furthermore, the proper time distribution of a vertex is independent of rank, which 
implies that the vertices tend to be aligned along a hyperbola in x — t space, cf. Fig ^ 

It is instructive (and historically more correct) to reverse the line of arguments. Searching 
for a left-right symmetric fragmentation function suitable for an iterative scheme, Eq (§) can 
be shown to be a unique solution. The final state density, G„, for an iteratively generated 
n-particle cluster with momentum fraction Z and mass is then 

(\.Z 1 z 

Gn{si, Z)dsidZ = 5(„(si)dsi — (1 - Z)" exp(-fei ^ ). (6) 

This density can be separated into one internal factor, independent of Z, and one external, 
depending only on Z and Si. The conventional choice is to let the internal part of the density 
be given by Qnisi), which then is to be interpreted as the final state density of a fragmenting 
string with mass ^/Sl. 

When we extend the model to several flavours, a set of new parameters enter. The inclusive 
probabilities Va that a quark has flavour a and the probabilities that a quark a and 
antiquark (3 form meson (or anti-meson) i are free parameters which satisfy the constraints 

a i 

It is also possible to have flavour specific parameters in the fragmentation function, which 
introduces a flavour dependence for the proper time distribution of a vertex. 

Assuming no flavour correlations between vertices, the left-symmetric fragmentation func- 
tion, Eq (H), generahzes to 

dz [ \ — 2 \ \)Tfl^ 

V(3G'^f,Mz,m^i)dz = vpG'^f,No,p{m^,)—z''- l^—j exp( -i), (8) 

which specifies the probability that flavour a is followed by flavour /3, forming a meson i 
which takes momentum fraction z. The N and a parameters are related by the unitarity 



3 



constraint J dzf{z,m'^) = 1. With the definition 



Cap{m ) = Nap{m )e 



2\ -bm^ 



ap\{bm 



(9) 



which in Eq ( |23| ) is shown to satisfy Cap = C^a, the density for an iteratively produced 
n-particle cluster starting with fiavour and ending with fiavour (3 is 



Cof3{Si) 



(10) 



The internal part of the density is given by 



r=lar,kr ^ 



P -bA 



'111 



which specifies the decay density of a string. 



The separation in internal and external parts is not unique, a factor depending only on si, 
tto and (3 can be associated with both parts. This introduces no arbitrariness to observables, 
depending on relative densities for different final states, as long as we adopt a factorization 
anzats, saying that each string, generated by parton-level physics, decays into exactly one 
hadronic final state. The density gli'^^ is always left-right symmetric in "truly" internal 
properties, i.e. quark fiavours Q;i...a„_i and hadron fiavours and momenta. In this paper the 
internal part is chosen to be symmetric also in its endpoint indices, g^'^^ = g^'"^ and to 
approach the simpler one-fiavour formulas in the limit of equal a^-values and equal masses 
rrii. 



We note that different phenomenological assumptions can be made to estimate the param- 
eters Va and G^^, and that, in general, several of them can be related to each other. In 
the standard Lund model, the string break-up is described by a tunneling process 0, which 
can be used to estimate suppression factors for strangeness and heavier spin states [|I|. In 
the UCLA approach the product vpGapNap is assumed to be proportional to spin- and 
isospin factors only, and the relative production of different mesons thus depend strongly 
on the factor exp{—bm'^/z) in the fragmentation function. (Strictly speaking, this implies 
flavour correlations between neighbouring vertices, i.e. the fa-parameters generalize to a set 
Vai3-) In this paper, there is no need to specify the flavour selection mechanism, and I will 
continue to use the general, though perhaps unnecessarfly numerous, parameters Va and G^^. 



3 MC and Model Differences 



Here I discuss implications of the differences in energy-momentum conservation between the 
Jetset MC and the Lund model, represented by Eq (p!l|). 
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Figure 2: Three final states. The densities for states II and III are determined by the same 
phase space factor, but state III is further suppressed by a large string area. However, with 
the standard MC energy conservation recipe, merging of a type I path from one string end 
and a type II path from the other can occur without additional suppression, thus enhancing 
the occurrence of type III states. 



In the iterative MC algorithm, hadrons are peeled off from the string ends until the remaining 
energy is getting small. Then two final hadrons are formed and their kinematics is essentially 
fixed by energy-momentum conservation. For 90 GeV strings, the approximate merging 
recipe affects two out of typically 10-15 primarily produced hadrons. For two-jet LEP events 
and experiments at more moderate energies, the final two particles constitutes a larger 
fraction of the primarily produced hadrons. 

The degrees of freedom indicated by the vague expressions "getting small" and "essentially 
fixed" are carefully specified in the MC to give the two merging hadrons reasonable prop- 
erties 0. Nevertheless, there are some differences. Consider the three final states in Fig ^ 
The high-multiplicity state I is enhanced by its phase-space factor but suppressed by a large 
area A. The low-multiplicity state II has a small area, but is suppressed by a small phase 
space factor. The state III is a permutation of II and, therefore, has the same phase space 
factor. Its larger area makes it significantly more suppressed than state //. 

Suppose that states / and // are equally probable. Letting both ends of the string fragment 
independently of each other implies that merging will occur between a /-type path from one 
side and a //-type path from the other just as often as between two paths of the same kind. 
This implies an enhanced occurrence of state ///. Thus, the fluctuations in vertex proper 
times is slightly larger in the MC, and the correlation between rank and rapidity somewhat 
reduced. 

In general, the MC and model differences are negligible. The bigger is the difference between 
state // and ///, the more suppressed are all the three discussed states. However, as 
mentioned in the introduction, this need not be the case at moderate energies well below 
the peak At high energies (like 90GeV) the energy conservation procedure may be 
of relevance for subtle enough observables. Within the string fragmentation picture, it is 
possible to make different assumptions about e.g. non-perturbative generation [|, |T^, |rT|, Q 



and baryon production [T^, ^ , giving different predictions for correlations between flavours, 
p± and rank. As rank is not an observable, a confrontation with data must exploit the 
correlation between rank and rapidity. We note that the merging procedure in Jetset, where 
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Figure 3: Production of n + m particles can be seen as production of two clusters. The 
hyperbolas in the figure show the possible end points for clusters with fixed mass squares 
Si and S2- There are two solutions, where the hyperbolas meet. Each solution has its own 
T-value, corresponding to the proper time squared of the vertex common to the clusters. 

the correlation between rank and rapidity is reduced 0, dilutes the observable differences 
between models for p±- and baryon generation. 



4 A new MC algorithm 

In this section I derive the new algorithm, working in the case with several flavours. After 
presenting the algorithm, I spend one sub-section on some technical details related to the 
computer time consumption. The section ends with a sub-section describing extensions 
needed to include transverse momentum generation in string fragmentation. 

We note that the decay density g{s) is zero or a 5- distribution for energies below the lowest 
two-particle threshold, cf. Eq (|1|). In the rest of this paper, I will constrain the discussion to 
energies above this threshold. 



4.1 1+1 dimensions 

We rewrite gn, Eq (|lTl), in the form 

i 
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(12) 



where 

A(a, b, c)=a^ + b^ + c^- 2ab - 2ac - 2bc, (13) 



r(a, 6, c) = - — 6 — c ± a/ A(a, 6, c) j . (15) 
Si, S2 and the two solutions to F are illustrated in Fig ^ 

The two-particle density g2 is directly calculable, as is also all meson-specific two-particle 
densities gij. For higher multiplicities we write 

9nS^) = 5^ t^,t;,G;,G^,^C,,(m,2)C,/3(m|) X 

dsids2gt''^\si) 



= J2 '"(^^ ^1, S2)^f'\s2, ml m]). (16) 
A(s,Si,S2)A(s2,m2,m2) u 

Eq (0) is the basis of the algorithm. Before describing it in detail, I comment on a few 
things. 



The cluster multiplicity is limited by n + 2 < vi/mmin- In general, the probability is 
negligible for much smaller n, since the typical multiplicity grows logarithmically with s. At 
this stage, it suffices to note that it exists an upper limit rimax to n, which implies that we 
can select a preliminary n from a fiat distribution. 

By ignoring the cluster momentum fraction Z, we can use iterative fragmentation to gener- 
ate the mass and endpoint fiavour of a n-particle cluster from the preliminary density (cf. 
Eq (0)) 

G^:'^\sMsi^ [ dZG("'^)(.i,Z)dsi = -^^(°'^)(si)dsi. (17) 

The denominator in Eq (|I6|) goes like [s2 — {rrii + rrij)'^]^ in the lower range of possible S2- 
values, and like [(v^— v^)^ ~ ^2]^ in the upper range. A preliminary distribution for S2, 
which takes these singularities into account and is suitable for MC generation, is 

i^„^ e(g2 - (m, + m,)^)e((yi-v/^)^ - ^2) ^^g^ 

TT ' [s2-(m, + m,-)2]^[(v^- 757)2-52]^ ' 
(s2 is easily picked from this distribution by setting 

S2 = ^ {(v^- v^)' + K + ^j? + [(v^- v^)' - K + ^j?] cos(i?7r)} , (19) 
where i? is a random number from a uniform distribution between and 1.) 
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After selecting n, Si and S2 as mentioned, the set of four possible final states obtained should 
be accepted with a probability proportional to the weight 

W^. ,5^$,$,, (20) 

[S2 - (mi - mj)2]2[(yi+y^)2 - S2]2 J J 

This weight has a finite maximum VFmax- The argument of C is in the finite range [m^j^, s], 
and there is thus and upper limit Cmax- This implies that 



34, 



3 

y / \ (ctmax ~'3'min) 



(21) 



W^max should be better estimated, to improve the acceptance of an event. That is discussed 
in the next subsection. At this stage, it suffices to note that it exists an upper limit Wmax, 
which implies that the acceptance probability is finite. 

The algorithm, including two-particle production, is as follows: 
1 Select n with the probability 
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n = 



(32+Wmaxn.max) ' 

Wmax 1 < r) < r) 

(n^-i-W n ^^> ^ — "' — "max 

\y2 max '^max j 



2a If n = select hadrons i and j with probability g^ij'^\s) / g^'^\s) and then solution I. 
After that, the fragmentation is completed. 

2b If n > iteratively generate n particles with the symmetric Lund fragmentation func- 
tion. Accept all momentum fractions Z for the generated cluster. This creates a cluster 
with the mass squared Si and end quark 7 from the density in Eq (|T7|). 



3 Select two final hadrons i and j with the probability "^^ ^7,5 ^5/3- < \fsi ^ 

rrii + rrij return to point 1. 



4 Select S2 from the preliminary distribution in Eq ([181). 

5 Accept the configuration with probability W/Wmax- Else, return to point 1. 

6 Select solution IJ with probability ^ j / ^^j^ j. Lorentz transform the cluster 
so that its last vertex is placed in the point specified by /. Set the momenta of the 
final two hadrons according to J . 



4.2 Computer Time Issues 

Above is presented an algorithm which, in finite time, produces a final state exactly according 
to the integrand in Eq (|l5). Many extensions are needed before the algorithm can give a 
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realistic description of hadronization. In particular, generalizing the algorithm to treat 
gluonic strings would require considerable effort. Before even considering such developments, 
it is relevant to discuss whether "finite time" also means "reasonable time" . In this subsection 
I estimate the time consumption with the new algorithm, as compared to an algorithm where 
no normalization numbers C need to be calculated, and where essentially each iteratively 
produced final state is accepted. 

The normalization numbers C must be numerically calculated, which in principle is a com- 
puter time problem. It is convenient to expand the integral representation of . The 
calculations are rather technical and not very illuminating, and here I only outline the pro- 
cedure. For details, cf. 0. 

The normalization number C for a cluster with mass ^/s, defined in Eq (P), can be written 
= y'||(M^"''-"'^^/^/'y^'^"(^)"^exp(-6.i^) = 
= , / da:d2/x''"i/"^exp(-x-i/)exp(-/). (23) 



Expanding exp(— |j) gives an asymptotic series which has to be truncated at some finite N, 
r-^l,] ~ ' V/ (« + a.,)!(« + af,)! 



As an example we get, for 0.5 and b = 0.75GeV"^ that this expansion can be 

used to obtain a precision 10~^ for all cluster masses above ~ 4.3GeV, and that it then 
suffices to include 13 terms in the expansion. 

For small bs, it is possible to expand (1 — z)°'i^ and do partial integrations to get an ex- 
pansion in bs. In the example 0.5 and b = 0.75GeV-^ it is enough to include 
approximately 70 terms in the bs expansion, to obtain a precision 10~^ in the whole range 
not covered by Eq (0). 

The sums in the bs expansion contain large numbers of opposite signs which almost com- 
pletely cancel each other. Thus, when very high precision is wanted, or when the a values are 
very large (significantly larger than the reasonable upper limit 1), higher precision than in 
the computer floating point types is required on each term of the sum. In such a situation, a 
special floating point representation, where the precision can be set by the user, is required, 
which implies longer computation times. However, this is not a problem in the example with 
reasonable a- values and precision demands discussed here, and the numerical calculation of 
C is not a big contributor to the computer time consumption. 

The efficiency E, i.e. the average probability to accept an event, is in the new algorithm 
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This number has been investigated in the simple case of only one hadron. The result is 
E ^ 1/30 for 10 GeV strings. As discussed in section H, there is little need to consider 
strings with higher energies. 

This efficiency holds when numbers nmax and VFmax, significantly better than the limits 
discussed after Eq (0), are found. At the start of a simulation, this is not the case, and 
computer time has to be spent on finding them. If identical strings are to be fragmented, 
the search for good rimax and VFmax need to be performed only once, but in general, when 
string configurations vary in different events, the calculation need to be done in every event. 
Results in the simple one-flavour case suggest that the flnal efficiency, including an empirical 
determination of rimax and PKnax in every event, is ~ 10~^. 

The conclusion of this subsection is that the new fragmentation algorithm is expected to be 
roughly a factor 100 slower than the standard Jetset MC. We note that this time compari- 
son applies only to the fragmentation phase. In a realistic event simulation, including parton 
cascades, resonance decays, detector simulation and event analysis routines, the relative time 
cost for the new algorithm will be substantially smaller. 



4.3 An algorithm with pj_ 

The model presented in section ^ needs to be extended to provide a realistic description of 
hadron production. Most extensions are briefly discussed in the outlook, but here I present 
in detail an algorithm including non-perturbative generation of transverse momenta. 

Adding transverse degrees of freedom implies that the two-particle density g2 no longer can 
be calculated exactly, but has to be generated by MC techniques. On the other hand, g2 is no 
longer singular, which means that two-particle states can be treated on more equal footing 
with the higher multiplicities. Thus, this model extension signiflcantly changes the algorithm, 
while many of the extra features discussed in the outlook should be more straightforward to 
implement. 

The standard way to describe k± generation in the string model is via a tunneling mecha- 
nism , which gives rise to an azimuthally symmetric Gaussian distribution. Here I present 
an algorithm applicable when a quark and antiquark in a vertex gets balancing transverse 
momenta k± from a distribution T{k±)d'^kj_ without singularities, and where there is no 
correlation between different vertices. 

Suppressing the k± dependence in the notation, the equation relevant for the algorithm reads 

[bs) 2 ^ J \'- ^ ) [ps2) 2 

where the density G^''^\si, Z)d?k±ry(lZdsi is obtained by iterative fragmentation. This will 
be used to determine the cluster variables si and Z, which implies that squared transverse 
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mass of the 2-particle state is 

S2 = {l-Z){s-si/Z). 



(27) 



Guided by Eq (p6D, we can formulate an algorithm as follows: 

1 Select n between and nmax from a uniform distribution. 

2 Generate iteratively n particles with the symmetric Lund fragmentation function. Re- 
turn to 1 if the remaining string energy is less than 2mmin. (if n = 0, no particles are 
generated. Instead set Z = 0, S2 = s and 7 = a.) It now remains to form two hadrons 
from the string remainder build up by a quark 7 with transverse momentum and 
an anti-quark (3 with transverse momentum —k±f^, and having a squared mass S2 given 



by Eq (13). 



3 Select flavour 6 for the final vertex and flavours i,j for the two final hadrons with 



probability vsGl^Gl 



4 Define perpendicular x and y axes so that ky^ — kyp = 0. Select kx in the last vertex 
from the distribution [/ dkyT(k±)~\ dk^- 



5 Calculate p^i = k^^ — k^ and p^j = k^ — k^p. If ^/s2 < a/^^ + + J w?- + ■ return 



to 1. 

6 Select a random number R from a fiat distribution between and 1. Set 



1 IX{s2,mf +pli,m] +pIj) 



ky = ky^ + ^d ^ ^ ^ 'J'' ' cos(i?7r) (28) 



7 Accept the event with a weight 



JdkyT{k^) {i-ZY^^2 ips,)-^ir[ ' ^ ^ ' 



If not accepted, return to 1. 
8 Select solution / and place the final two particles accordingly. 



5 The Total String Decay Density 

As discussed in section ^ observables depend on relative densities for different hadronic 
final states possible from a string, and not on its total decay density g{s). Nevertheless, 
a numerical calculation of g{s) is of interest. It is instructive to see at what energies g{s) 
— in whatever form it has been specified — approaches its asymptotic behaviour. When 
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Figure 4: The decay density g{s)/b as a function oj ^fsjm, jorhm^ = 3/8 and three different 
a values. At = 2,3,4 clear threshold effects are seen, but already at ^/s ~ 6m the 

density has well approached its asymptotic form, indicated by dotted straight lines. Results 
for other values of bm? , corresponding to meson masses in the range 0.5 to 1 GeV, are 
similar. 



the asymptotic limit is valid, the iterative algorithm, postponing the energy momentum 
conservation problem to a later stage, is a good approximation. 

It is much simpler to calculate the densities gn{s), and thus their sum g{s), than to generate 
final states. In the 1 + 1 dimensional model without p± generation it can be done as follows: 
The two-particle density, g2{s), is directly calculable. For all higher multiplicities n + 2, 
generate an n-particle cluster and calculate its mass ^/sl. Set the "merging weight" W = 
if it is kinematically impossible to produce a final state with a cluster of mass y/si and two 
more hadrons. Else, select a mass y/s2 for the final two particles from the distribution in 
Eq (^) and calculate the weight W in Eq (|2Dp. Repeating the procedure gives gn{s) as the 
average of W. 

Fig ^ shows results in the simple model without p± and with only one meson fiavour of mass 
m. The dimensionless quantity g{s)/b is plotted as a function of y/s/m, for different values 
of a. The quantity bm"^ is set to 3/8, corresponding to a typical value 0.75GeV~^ on b and 
typical mass m = 0.7GeV. One can clearly see the threshold effects for = 2, 3, 4. After 

that, the contribution of gn{s = n^rn^) is getting less important relative to lower multiplicity 
contributions, and the sum g{s) is well approximated by its asymptotic form at ^/s ~ 6m. 
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This implies that high energy strings may be iteratively fragmented using the symmetric 
fragmentation function, as long as there is a negligible risk to get a string remainder with 
an energy below some six hadron masses. In doing so, the algorithm efficiency, discussed in 
section |4.2| for a string energy of about 15m, does not deteriorate with increasing energy. 



6 Summary and Outlook 

The Lund string fragmentation picture [jl| is successful in comparison with data, and the 
computer implementation Jetset is a widely used event generator for high energy physics 
experiments. However, there are differences between Lund String fragmentation model and 
Monte Carlo. Global energy momentum conservation does not fit into an iterative fragmen- 
tation scheme and is observed in a slightly different way in the MC, as compared to the 
model mathematical expressions. The differences may be non-negligible in certain cases, like 
few-body production p and subtle enough observables, as those discussed in [Q. 

I present an algorithm which observes energy-momentum conservation exactly as proposed 
by the model expressions. An embryo to a computer program based on this algorithm ^ is 
used to calculate the total decay density g{s) of a string. It is shown to be well approximated 
by its asymptotic behaviour for string energies above some six hadron masses. 

A complete fragmentation simulator has not yet been developed, but only the parts needed to 
calculate g{s). With this program embryo, it is possible to estimate the time consumption of a 
future event simulation program, which is found to be roughly 100 times larger than Jetset. 
It should be noted that this is only in the fragmentation phase. In a realistic calculation, 
including parton cascades, resonance decays, detector simulation and event analysis routines, 
the relative time cost for the new algorithm will be substantially smaller. 

A set of extensions to the algorithm are needed, before a realistic generation of hadronic final 
states are possible. Inclusion of fiavour- and p± correlations as in |]TU|, baryon production 
as in 



jT2[, modifications to the area-definition in the case of heavy endpoint quarks 
and Breit-Wigner mass distributions for unstable resonances should all be straightforward 
developments, while more advanced baryon models and strings with gluons may need 



more effort. The latter is a particularly complex problem, as apparent from the thorough 
description in [|T5|. 
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